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Abstract The present paper reviews recent advances in the theory of nonlinear driven 
magnetohydrodynamic (MHD) waves in slow and Alfven resonant layers. Simple es- 
timations show that in the vicinity of resonant positions the amplitude of variables 
can grow over the threshold where linear descriptions are valid. Using the method of 
matched asymptotic expansions, governing equations of dynamics inside the dissipa- 
tive layer and jump conditions across the dissipative layers are derived. These relations 
are essential when studying the efficiency of resonant absorption. Nonlinearity in dis- 
sipative layers can generate new effects, such as mean flows, which can have serious 
implications on the stability and efficiency of the resonance. 
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1 Introduction 

The dynamical response of the plasma in the solar atmosphere to rapid changes can 
be manifested through wave propagation along and across the magnetic field. Many of 
these waves are in the magnetohydrodynamic (MHD) threshold (waves with periods 
larger than the ion collisional time and wavelengths larger than the mean free path 
of ions) and they are observable by the new generation of space- and ground-based 
telescopes in almost all regions of the solar atmosphere (for extensive reviews on wave 
obser vation see, e.g.. lActon at al.1ll98ll ; lNakariakov and Verwichtel2005l : lBaneriee et al.l 
120071 ). 

One of the most fundamental characteristics of solar and space plasmas is their 
very high degree of inhomogeneity along and across magnetic fields. It is well known 
that the properties of MHD waves are strongly modified by the plasma inhomogeneity. 
In particular, when the inhomogeneity of the medium is transversal to the direction of 
wave propagation a new phenomenon, called resonant absorption, can appear. 
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While homogenous plasmas have a spectrum of linear eigenmodes which can be 
divided into slow, fast and Alfven subspectra, with the slow and fast subspectra hav- 
ing discrete eigenmodes and the subspectrum of Alfven waves being infinitely degen- 
erated, in inhomogeneous plasmas the three subspectra are substantially changed. 
Acc ording to spectral theories of waves in inhom ogeneous plasmas (see, 
e.g.. lGoedbloedlE98i : lGoedbloed and Poedtslfiooih the spectra of slow and 
Alfven waves become continuous while the spectrum corresponding to fast 
waves becomes discrete. Eigenfunctions that correspond to frequencies in the con- 
tinuum spectra are improper as they contain a non-square integrable singularity at the 
resonant position. 

According to the accepted resonant wave theories, effective energy transfer be- 
tween an energy carrying wave and the plasma occurs if the frequency of the wave falls 
into the slow or Alfven continuum, i.e. at the slow or Alfven resonances. 

The process of resonant interaction of waves and the possibility of energy transfer 
between global and local waves has been studied intensively for the past few decades. 
First resonant absorption was studied as a means of supplemen tary heating fusion plas- 
mas, but was later rejected due to technical difficulties (see, e.g., Tataronis and Grossmanl 



1 19731 ; IChen and Hasegawdfl97l IPoedts et aT1ll989l ; Ivan Eester et al.lll99ll V 



In the Earth's magnetosphere resonant MHD wave coupl ing is believed to gen- 
erate low frequency pulsation or energize ULF waves (see, e.g. Lanzerotti et al.l 19731; 



Southwoodll974l;ISouthwood and Hughesjl983l:lRuderman and Wrightll99*^ Erdelvi and Tarovanl 



2003: iTarovan and E rdelvi 2002, 2003a|(bO. Within the context of solar physics Ionson 



(1978) was the first who proposed resonant absorption as a possible mechanism for 



authors (see. e.g. iKuDerus et all 


1983: [Bavilal 19871: Hollweg 199ll: Goossensl 


1991; 


Ofman and Davilalll995l: Erdelvi and Goossensll 19941. Il995. 19961: Erdelvill 19971 


1998; 


Ballai et al. |2000|. etc.), however it became clear soon that resonant absorption alone 



cannot explain the very high temperature of the corona. Resonant absorption was 
also used to explain the loss of power in p-modes when i nteracting with sunspots in 
the solar photosphere (see, e.g., Ilfollweg 1988|: Loul ll990l; Goossens and Poedts Il992l; 



Spruit and Bodganl[l992llErdelvi and Goossensj|l994l ; lKeppens et al.ll994lErdelvill9'96l ; 



Tirrv et al. lll998d ; lTfrrv et alj|l998bh . 

Recently resonant absorption acquired a new application in coronal seismology 
where the rapid damping (with damping times of the order of a few periods) of kink 
oscillations in coronal loops is expl ained in terms of resonant interaction of global kink 
modes and the local Alfven waves jRuderman and Robertsll2002l ; iGoossens et alj|2002 

i 
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Dvi 



Terradas et al 



and Rudermanll2006l ; lErdelvi and Verthll2007l ; IVerth and Erdelvill20Q 
20101 ). Coronal seismology is dependent on theoretical relations (dis- 



persion relations) which link plasma parameters, such as the plasma density, to wave 
parameters, such as the wave frequency, in a precise way. Generally, plasma parameters 
are determined from wave parameters, which themselves are determined observation- 
ally. The dispersion relations for many plasma structures under the assumptions of ideal 
MHD are well known; they were derived long before accurat e EUV observations were 
available fe.g. lEdwin and Roberts! 1 19831 ; [Roberts et al.lfl984T) using simplified models 
within the framework of ideal and linear MHD. Although the realistic interpretation of 
many observations are made difficult by the spatial and temporal resolution of present 
satellites not being adequate, considerable amount of information about the state of 
the plasma and the structure and magnitude of the coronal magnetic field have been 
already obtained. 
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At resonance, amplitudes of oscillations (i.e. the energy density) can grow without 
limits, however even a small amount of dissipation is able to prevent the unlimited 
increase of the amplitude. The presence of dissipation causes transforming the wave 
energy into heat. On the other hand solar and space plasmas are media with very high 
Reynolds numbers (i.e. weakly dissipative media), which means that the dissipation 
is an important ingredient in the description of dynamics only in a narrow region 
around the resonant magnetic surfaces called dissipative layer. Outside this layer the 
dynamics is described by the ideal MHD. This natural structuring of t he domain of 
interest allows us to use th e method of matched asymptotic expansions (|Navfehlll98ll : 
iBender and Orszaghl [l99lh to describe the problem of resonant absorption of MHD 
waves. The essence of this method resides in solving ideal MHD equations outside the 
dissipative layer, the non-ideal equations inside the dissipative layer (which contains 
the resonant surface), and matching the two sets of solutions in overlapping regions 
(for details see the review by Goossens et al. 2010, this issue). 

The dissipation of energies in the solar atmosphere is a delicate problem as the 
dominance of a certain process over all other possible ones depends on the region 
where the dynamics will be described and also on the very physics of the problem 
under consideration. That is why the dynamics of waves in the solar atmosphere is 
going to be influenced by different dissipation mechanisms in, e.g. the photosphere and 
solar corona. A full analysis of particular dissipative mechanisms used in the present 
paper will be given in the next section. 

Nonlinearity was systematically overlooked in many studies related to wave dy- 
namics because the mathematical treatment of nonlinear problems is cumbersome, 
and because the complicated mathematical analysis can often obscure the physics and, 
therefore, make the results unusable for observers. Dynamics occurring over extended 
scales (as in the solar atmosphere) are likely to encounter situations which can lead to 
the increase in amplitudes resulting in breakdown of linear description. Very often the 
wave amplitudes can increase due to the change of the medium properties, e.g. acoustic 
waves propagating from the photosphere to the corona can steepen into shocks due to 
the decrease of the density with height. 

A great advantage in tackling the complicated problem of r esonant absorption i s 
the notion of connection formulae introduced for the first ti me by Sakurai et alj ( 199ll ) . 
The accuracy of the method by numerics was addressed bv lStenuit et al .1 <|l995h . This 
approach assumes that the thin dissipative layer is treated as a surface of discontinuity 
with the dynamics at both sides of the discontinuity fully described by the ideal MHD. 
The system of non-ideal MHD equations describing the dynamics inside the dissipative 
layer is used to obtain the connection formulae that serve as boundary conditions at 
the surface of discontinuity in the same way as the Rankine-Hugoniot jump conditions 
for shocks. 

A great deal of understanding of the process of r esonant absorption 



came from recent numerical i nvestigations (see, e.g.. lOfman at "all Il994l : 
lOfman and Davilalll995l . 1 19961 : IPoedts and Goedbloedl Il997h . where issues 
like higher dimension resonance, time evolution of the resonance or ran- 
domly driven resonance were discussed in great detail. 

The paper is organized as follows. In Section. [2] we introduce the equations used 
to describe the dynamics of waves and present the mathematical and physical tools 
employed in the mathematical presentation of the problem. In Section 3 we introduce 
scalings and dimensionless quantities required to quantify the importance of nonlin- 
earity. The problem of nonlinear slow resonant waves is presented in Section 4, where 
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we describe the nonlinear character of waves near resonance and we provide the jump 
conditions near the resonance. In Section 5 we show that there are two particular cases 
where the jump conditions can be simplified into an explicit form. The problem of 
nonlinear resonant Alfven waves is discussed in Section 6. As an application to the 
nonlinear theory reviewed in this paper we will apply the equations to study the ef- 
ficiency of resonance in the limiting cases of weak and strong nonlinearity. A direct 
consequence of the nonlinear framework used in our analysis is the generation of mean 
shear flows at both resonances, which will be presented in Section 8. 



2 Basic equations 

Solar and space plasmas are far from being an ideal environment with the plasma 
dynamics being affected by many dissipative and dispersive effects. As stressed earlier, 
the right choice of dissipative mechanism depends on the location where a particular 
physical phenomenon takes place as well as on the nature of this phenomenon. For 
the present paper, we are going to consider viscosity, electrical resistivity and thermal 
conduction. Dissipative processes are weak in the solar atmosphere; this means that 
the diffusion coefficients are small. The rate of dissipation, however, is dependent on 
the local spatial scales. 

Many phenomena (e.g. resonant absorption, current sheets or turbulences) are in- 
herently non-ideal and very often nonlinear as they are strongly influenced by dissi- 
pative and/or dispersive effects. In particular, dissipation is important for nonlinear 
dynamical processes because large-scale motions rapidly lead to the formation of small- 
scale structures, which corresponds to singularities in the ideal theory. 

Viscosity and thermal conductivity are linked to hydrodynamical processes, while 
electrical conductivity and Hall dispersion are related to the presence of the electric 
and magnetic fields. In the context of solar physics, the general effect of dissipation and 
dispersion is to relax the accumulation of wave energy in a system. The relaxation can 
be performed by, e.g. converting the wave energy at resonance into heat by viscous or 
resistive dissipation or thermal conduction (via resonant absorption) , or the dispersion 
of energy over a larger area by the Hall effect. The relaxation caused by dissipation 
and dispersion prevents the formation of singularities (entities abhorred by nature). 

The key quantities in our discussion are the dimensionless products o;,;t; and cj e T e , 
where <^i( e ) is the ion (electron) cyclotron frequency and T~i( e ) is the mean ion (electron) 
collision time. Since viscosity is mainly due to ions, the product WjT^ is important in 
describing viscosity. On the other hand, thermal and electrical conductivity are effects 
attributed mainly to electrons, so the product ui e T e appears in the description of these 
two dissipative mechanisms . 

The classical iBraginskiil (|l965h expression for the viscosity tensor reads 



S = £) W S i , (1) 

i=0 

where the coefficients 770, r\\ and 772 describe viscous damping and 773 and 774 are non- 
dissipative coefficients related to wave dispersion due to the finite ion gyroradius. The 
quantities Sj are given in terms of the unit vector along the equilibrium magnetic field 
and the velocity vector. The viscous force in plasmas is given by V • S. In general the 
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expressions for % (i = 1, . . . , 4) and Si are quite complicated and we do not write them 
down. The coefficient of compressional viscosity r/Q is given by 

pk B TTj 

Vo ~ 1 2 

m p 

where p and T are the density and temperature of the plasma, kg is the Boltzmann 
constant and m p is the proton mass. When oj,;t.; <C 1 as in the solar photosphere, the 
other coefficients in Eq. {1} are given by 

m-m-vo, m ~ 2774 « 0.8770(0;,^). (3) 

It can be seen that 773 , 774 <C 770 , so that the last two terms in Eq. {TJ can be neglected. 
Then the viscous force is given by the approximate expression 

F vis = V-S^77o(V 2 v+|vV-v), (4) 

where v is the plasma velocity. This expression for the viscous force is isotropic in the 
sense that the force is independent of the magnetic field direction. 

In the solar chromosphere and, especially in the corona, the condition cjjTj 3> 1 is 
satisfied. In this case the remaining four dissipative coefficients are given by 

^ = 77~^2 ' ^ 2 = 4r ?l' ? ?3 = m = 2?73- (5) 

Now 771,772,7/3,774 <C 770, so that, in particular, the compressional viscosity strongly 
dominates the shear viscosity. As a result all terms in Eq. |l} can be neglected in 
comparison with the first one, and we arrive at the following approximate expression 
for the viscous force, 

F vis = V ■ S w 770V ■ (b ® b - §f) [3b ■ (b • V)v - V • v]. (6) 

Here bo is the unit vector in the magnetic field direction, I is the unit tensor, and f ® g 
denotes the dyadic product of vectors f and g. 

However, Eq. © can be used only in slow dissipative layers. The compressional 
viscosity is identica lly zero for Alfven waves. As a result, i t cannot remove the ideal 
Alfven singularity (lOfman at all 1 19941 : lErdelvi et al.l 1 19951 : iRuderman and Goossensl 
Il996l :l Mocanu et al.l I2OO8I ) . When studying the motion in Alfven dissipative layers we 



have to keep the second and third terms in Eq. (JTJ). On the other hand, in spite that 
Wii i]l <K 770, we can neglect the first term in Eq. ([1]). As for the fourth and fifth terms, 
they can be neglected because they do not provide dissipation. In addition, we can 
neglect the spatial variation of the unit vector in the magnetic field direction, b, in the 
Alfven dissipative layer. As a result, we obtain the following expression for the viscous 
force, 

F vis w 771V • (Si + 4S 2 ) w 771 {V 2 v + 3b V 2 (b ■ v) + 4b b ■ VV ■ v 
+ 3(b ■ V)[(b ■ 3V) + V(b ■ v)] + V[b • V(b ■ v)] 
- 15b b -V[(b -V(b -v)]}, (7) 

which is still rather complicated. However, in Alfven dissipative layers V • v « which 
enables us to neglect the third term on the right-hand side of Eq. (J7J. In addition, in 
Alfven dissipative layers large gradients can develop only in the direction perpendic- 
ular to the equilibrium magnetic field. This observation enables us to neglect also the 



6 



fourth, fifth and sixth terms on the right-hand side of Eq. (01 . Finally, the dominant 
component of the velocity in Alfven dissipative layers is the component orthogonal to 
the equilibrium magnetic field, while the component along the magnetic field remains 
small. This means that we can also neglect the second term on the right-hand side of 
Eq. @. As a result the very simple expression for the viscous force that can be used 
in Alfven dissipative layers is 

F vis » m V 2 v. (8) 

The wave energy can be also dissipated due to thermal conduction. In what follows 
we assume that the ions end electrons have equal temperatures. In this case the general 
expression for the heat flux q is given by 

q = -K|| VjjT - k_lV_lT - K A b x VT, (9) 

where K||, k±_ and ka are the parallel, perpendicular and skew coefficients of thermal 
conduction, V|| = bo (bo ■ V), Vj_ = V — bo (bo • V), and T is the plasma temperature 
related to the plasma density, p, and pressure, p, by the ideal gas law, 

P = ^ P T. (10) 

For the electron-proton plasmas, the expression for km is given by 



3kgpTr e 
rripm e ' 



(11) 



where m e is the electron mass. When uj e T e < 1 as in the lower part of the solar pho- 
tosphere, kj_ w km and ka ~ (cj e T e )K| <C kii. As a result the approximate expression 
for the heat flux reads 

q=-K||VT. (12) 

This expression is isotropic in the sense that it is independent of the magnetic field 
direction. 

When cu e T e 2> 1 as in the solar chromosphere and corona, ~ (^ei~e) 2 K|| <IC km, 
ka ~ (o; e r e )~ 1 K|| <C km, and the first term in Eq. @ strongly dominates over two 
other terms. Therefore we can use the approximate expression for the heat flux 

q=-K||V||T. (13) 

We now see that the heat flux is in the magnetic field direction, and magnetic surfaces 
act as thermal insulators. 

We write the expression of the generalized Ohm's equation in the form 

E + vxB = - + ^jxB, (14) 
a pe 

where e is the elementary electrical charge, a the electrical conductivity, j the electrical 
current density, and E and B are the electric and magnetic field, respectively. The 
first term on the right-hand side of Eq. (|14p is related to the plasma resistivity. It 
disappears in the limit of infinitely conducting plasmas where a — > oo. The second term 
on the right-hand side of Eq. (|14|) describes the Hall current related to the account 
of the ion inertia. A more general form of the Omh's equation contains also the term 
proportional to B x (j x B) related to the so-called Cowling conductivity (see, e.g. 
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Priest 2000). This term is identically zero for a fully ionized plasma, while in partially 
ionized plasmas, can be important only when the electron gyrofrequency is large and 
the plasma is sufficiently rarified, so that the mean collision time between electrons and 
neutrals is large. On this ground we disregard this term. In what follows we neglect 
the displacement current in Maxwell's equations and use Ampere's law to relate the 
electrical current density and magnetic field, 

M j = V x B, (15) 

where no is permeability of free space. 

In our analysis we use the system of MHD equations 

g+V-pv = 0, (16a) 



p (§5T + v ' Vv ) = _Vp + x B) x B + Fvis ' (16b) 

f-VxE, (16c) 

In these equations F v ; s is given either by Eq. @, Eq. (O or Eq. ((SJ), and E is given by 
Eq. (|14[1 ; 7 is the ratio of specific heats or adiabatic index, and C is the loss function 
given by 

j 2 

£ = V- v- =--S:Vv, (17) 
a 

where S is the viscosity tensor and the colon indicates the double summation, 

S:Vv S,,±-L. (18) 

Recall that the magnetic field is solenoidal, i.e. 

V-B = 0. (19) 

This equation can be considered as an initial condition imposed on the magnetic field 
because, if it is satisfied at t = 0, then it follows from Eq. (|16cp that it is satisfied for 
any t > 0. If we neglect the Hall term in Eq. (|14|l then the induction equation (|16c|l 
reduces to 

^ = Vx(vxB) + nV 2 B, (20) 
at 

where 77 = I/hqct is the coefficient of magnetic diffusion. These equations must be 
supplemented by the Ohm's law (I14[) expressing the connection between 
the electric field and magnetic field. 
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3 General discussion and dimensionless parameters 

Traditionally the importance of dissipative processes in MHD is characterized by di- 
mensionless parameters calculated as ratio of corresponding dissipative terms to dy- 
namical terms. The importance of viscosity is quantified by the viscous Reynolds num- 
ber. In the previous section we have seen that the coefficients of compressional and 
shear viscosity can be quite different, so that we have to introduce two Reynolds num- 
ber, one related to compressional and another to shear viscosity, 

R% = !±™, R% = !±™ t (21) 

■no m 

where V is a characteristic speed, usually taken to be equal to the phase speed of waves, 
L a characteristic length, and p* a characteristic density. Of course, where the viscosity 
is isotropic, R% = R% = R e and F v i s is given by Eq. Q. The characteristic length L 
plays a very important role in Eq. (|21 [) . Far from resonant layers it can be taken to 
be equal to the wave length. Typically in the solar corona we obtain very large values 
of R% and Rf., which enables us to neglect viscosity, for example, R% = I0 3 -j- 10 4 and 

R% = I0 12 -r I0 14 . 

The importance of resistivity is characterised by the magnetic Reynolds number 
defined as 

Rm = — . (22) 

V 

For the solar corona typically R m ~ 10 12 -j- I0 14 . 

Finally, the importance of thermal conduction is characterised by the Peclet number 
which is given by 

p e = hBPiXL. (23) 

m p K|| 

This number is very large in the solar photosphere (I0 8 — I0 9 ), but it can be of the 
order of 100 in the solar corona. 

The situation with the resistivity and thermal conduction is similar to that of 
viscosity. Usually we can neglect them far from dissipative layers and consider plasmas 
as infinitely conducting (p — > oo) and their motions as adiabatic. However, these two 
dissipative processes can be again very important inside the dissipative layer. 

Linear theory p redicts the existence of slow and Alfven dissipative layers (e.g. 
iGoossens et aljfeoioh . When the dominant dissipative processes are the isotropic plasma 
viscosity and resistivity, as in the solar photosphere, it is convenient to introduce the 
total isotropic Reynolds number, Ri, defined by 

Ri Re Rm ^ 

Let us also introduce the dimensionless wave amplitude far from a dissipative layer, 
e -C 1. Then linear theory predicts that the amplitudes of perturbations of certain 
variables called 'large variables' in a slow dissipative layer embracing the slow resonant 

1/3 

position are of the order of tR- . The linear theory can be used to describe the motion 
in the slow dissipative layer only when dissipation in this layer strongly dominates over 
nonlinearity. Estimates show that the ratio of the largest typical nonlinear term in the 
MHD equations to the largest dissipative term is of the order of (see, e.g. Ruderman 
et al. 1997a) 

N i = eR 2 i /3 . (25) 
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The linear description is only valid when JV, -C 1. In the solar photosphere, however, 
6 = 1CT 2 -j- 1CT 3 and i?j > 10 6 , so that JVj > 10. This estimate clearly shows that the 
motion in slow dissipative layers is strongly nonlinear in the solar photosphere. 

In the solar corona the dominant dissipative processes affecting slow waves 
are the compressional viscosity and thermal conduction (recall that it is strongly 
anisotropic and the heat flux is parallel to the magnetic field). To characterise both dis- 
sipative processes simultaneously we introduce the total anisotropic Reynolds number 
defined by (see, e.g. Ballai et al. 1998a) 

7T = 7^ + F- (26) 

tla ri e r e 

The linear theory predicts that the amplitude of large variables in dissipative layers is 
of the order of eR a , while the ratio of the largest nonlinear term in the MHD equations 
to the largest dissipative term is of the order of 

N a = eR 2 a . (27) 

Again, the linear description is only valid when N a < 1. In the solar corona for a 
typical e = 10" 2 10~ 3 and R a > 100, so we obtain N a ^ 10, and the motion in slow 
dissipative layers in the solar corona is strongly nonlinear. 

Viscosity and resistivity are also the dominant dissipative processes in Alfven dis- 
sipative layer in the solar photosphere. Hence their importance can be characterized 
by Ri. Since the quasi- Alfvenic motion in Alfven dissipative layers practically do not 
perturb the plasma temperature, thermal conduction is not important in these layers. 
Since, in addition, compressional viscosity also does not operate in Alfven dissipa- 
tive layers, the dominant dissipative processes in these layers in plasmas with strongly 
anisotropic viscosity like in the solar corona are the shear viscosity and resistivity. 
These two processes can be simultaneously described by the total Alfven Reynolds 
number Ra defined by an equation similar to Eq. (|24[) but with R% substituted for R e . 
Applying the same formal procedure as that resulted in Eq. (|25[1 we obtain a similar 
estimate for the ratio of the typical largest nonlinear term in the MHD equations to the 
largest dissipative term N^. On the basis of this estimate we conclude that the linear 
description of Alfven dissipative layers breaks down as soon as e > R^^ . However, 
as it will be explain later, this conclusion is wrong as in Alfven dissipative layers the 
largest nonlinear terms in the MHD equations cancel each other. 



4 Governing equation for motion in slow dissipative layers 

In what follows we adopt Cartesian coordinates x, y, z and assume that the equilibrium 
quantities depend on x only, while the equilibrium magnetic field is unidirectional and 
parallel to the j/z-plane. In the equilibrium plasma is at rest. We only consider a two- 
dimensional problem and assume that perturbations of all quantities are independent 
of y. 

In order to derive the equation governing the wave motion inside slow dissipative 
layers we need to use the matching conditions with the solution of linear ideal MHD 
equations outside the dissipative layer, which we will call the external solution. Hence, 
before we embark on the derivation of the governing equation for slow dissipative layer 
we briefly discuss this external solutions. 
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To obtain the external solution we write all variables in the system of MHD equa- 
tions as a sum of an equilibrium quantity and their Eulerian perturbation 

p = p +p', p = p Q + p', B = B + b. (28) 

We substitute these expressions in the ideal MHD equations and linearize them with 
respect to perturbations. We Fourier-analyse the obtained linear equations with respect 
to z and t and take all quantities proportional to exp(ifcz — itot). Finally, we eliminate all 
the variables from these equations in favour of the s-component of velocity, u, and the 
total pressure perturbation P = p + Bo ■ b//j,o to obtain the system of two first-order 
ordinary differential equations of the form 



du iuiP dP ikpoDj^u 
dx^~F~' ~dx = V ' 

In these equations V = ui/k is the phase velocity and 



(29) 



V 4 — V 2 (v A + c|) + v\<? s cos 2 a V 1 

Da = V 2 — v\ cos 2 a, Dc = (i>a + c s)(^ 2 ~~ C T cos2 Q )i (31) 
where the Alfven va, the sound c$, and cusp ct, speeds are defined as 

2 _ 2 _ IPO 2 _ v 2 a c S 

V A - 1 C S - 1 C T - — 2"' ( 62 > 

A'OPO PQ v A + Ct, 

with a being the angle between Bo and the z-axis, i.e. cos a = Bq z /Bq. 

The system of equations (|29l) nas two regular singular points, xa and x c . The 
first point is defined by the equation Da(xa) = and corresponds to the Alfven 
resonance. The second point is defined by the equation D(j(x c ) = and corresponds 
to the slow or cusp resonance. Using the Frobenius expansions one can show that the 
asymptotic behaviour of solution in the vicinity of these points is given by P ~ const 
and u ~ const x In \x — xa c |. Hence, Pis a regular function in the vicinity of the Alfven 
and slow singularity. The asymptotic behaviour of other quantities can be determined 
from the linearized MHD equations. Hence, in the vicinity of the Alfven resonance, 
this behaviour is given by (see, e.g. Sakurai et al. 1991; Ofman et al. 1994; Erdelyi and 
Goossens 1995; Erdelyi 1997, 1998) 

const 





j const x In \x - 


■ x A \, 


v \\ 


~ const, 


v± 


b x 


~ const x In \x ■ 




b r 


~ const, 


bx 


,>'■ 


••> const, 




p>, 


^ const, 


P ■ 



x — xa 
„ const ; (33) 
x — xa ' 
const, 



where v» is the velocity component parallel to the equilibrium magnetic field, and vj_ 
is the velocity component perpendicular both to the equilibrium magnetic field and to 
the ^-direction. These two quantities are defined by 



'II 



Bo/Bo, «l = v ■ (B x e x )/B , (34) 
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where is the unit vector in the ^-direction. The quantities 6|| and bj_ can be deter- 
mined in a similar way. The asymptotic behaviour of perturbations at slow resonance 
are (see, e.g. Sakurai et al. 1991; Ballai et al. 2000a) 

u ~ const x In \x — x c 



const 






i 


V± 


~ const 


X Xq 




const 


b± 




i 


~-> const 


X Xq 




const 






i 


P - 


const. 


X — Xc 







I) "II 

b x ~ const x In \x — x c \, bu ~ , b± ~ const, (35) 

/ const / 



X Xc 

The obvious result of Eq. (|33[1 is that the most singular variables at x = x& are v±_ 
and bj_ and here they have a l/i singularity. Due to their behaviour at the resonant 
position they are called the large variables. In accordance with Eq. (|35|1 . the large 
variables at the slow resonant position x c are vu, bu, p' and p' . 

An important consequence of Eqs. (|33|) and (|35|l is that the perturbation of the 
total pressure, P, does not change across the dissipative layer. The physical reason 
for this behaviour resides in the fact that since the dissipative layer is thin, i.e. its 
inertia is very small, the acceleration of the dissipative layer is finite, the total force 
applied to it has to be small. This implies that the total pressure at the two sides of the 
dissipative layer has to be almost the same. If we denote the thickness of the dissipative 
layer as t and the characteristic length of the problem as L, then the variation of the 
total pressure across the dissipative layer is zero in the zero-order approximation with 
respect to the small parameter £/L. The conclusion that there is no jump of total 
pressure perturbation across the dissipative layer is obtained here on the basis of ideal 
MHD. This conclusion is confirmed by the linear theory of dissipative layers based on 
dissipative MHD (e.g. Goossens et al. 1995 for Alfven resonance and Erdelyi 1997 for 
slow resonance). 

The linear theory suggests the following physical picture of wave motion in slow 
dissipative layers. The global motion of the plasma is in exact resonance with slow 
waves at the ideal resonant position, and it is in quasi-resonance with slow waves in a 
narrow dissipative layer embracing the resonant position. The variation of the external 
total pressure acts as a driver of slow waves in the dissipative layer. Hence, the aim 
of nonlinear theory is to derive the equation governing the nonlinear evolution of slow 
wa ves in the dissipative laye r. 

iRuderman et al] ( 1997ah were the first who considered this problem. In their anal- 
ysis they assumed that the only dissipative processes operating in the plasma are resis- 
tivity and isotropic viscosity, and both dissipative processes are weak, i.e. dissipation 
is characterized by Ri S> 1. They considered the motion in the form of a propagating 
wave with permanent shape, so that perturbations of all variables depend on the linear 
combination 9 — z — Vt rather than z and t separately. The phase speed of slow waves 
at x c is equal cx c cos a, where the subscript V indicates th at a quantity is calculated at 
x — x c , so V — cxc cos a. Finally. [Ruderman et alj ( 1997al ) assumed that perturbations 
of all va riables are periodic with respe ct to with the period (wave length) equal to 
L. Later IRuderman and Erdelvil (|2000h generalized this derivation to allow slow time 
variation of the wave shape, so that perturbations of all variables are functions of x, 9 
and 'slow' time r. In what follows we briefly outline this derivation. 

To obtain the governing equation that simultaneously describes both nonlinear and 
dissipative effects we formally take Ni ~ 1, so that Ri ~ e~ 3 / 2 . This assumption is 
formal in the sense that one can consider iVj <C 1, in which case we can neglect non- 
linearity in comparison with dissipation, and iVj <C 1, in which case we can neglect 
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dissipation in comparison with nonlinearity. Linear theory predicts that the character- 
istic thickness of the dissipative layer is LR^ S ~ e 1 / 2 . This implies that it is convenient 
to introduce the stretching variable £ = e~ 1 ' 2 (x — x c ) in the dissipative layer. It is also 

1 /n 

practical to introduce the 'slow' time r = e ' t. 

The nonlinear interaction of the wave motion in the dissipative layer with the 
plasma causes the mean flow in the yz-plane with the amplitude of the order of e ' . 
Hence, it is convenient to split the velocity components in the mean and oscillating 
parts, 

V\\= «||+«||» f_L='"±+£_L, «|| = («||>) «L = («-L>> (36) 
where the mean value of a period function /(#) over the period is defined by 

It follows from Eq. (|35|) that, in the dissipative layer, vu ~ e 1 ^ 2 , 6m ~ e 1 ^ 2 , p' ~ e 1 ^ 2 , 
p ~ e 1 ^ 2 , u ~ elne, frz ~ elne, ~ e, ~ e, and P ~ e. Since |lne| is much 
smaller than any negative power of e for e <C 1, in what follows we consider me as 
a quantity of the order of unity and we look for the solution to the dissipative MHD 

1 /2 

equations in the form of power series expansions with respect to e ' . We write this 
expansions in the form / = e^/i + e/2 + • • • for p' , p , vu , 61 , vu and vj_, and in the 

form g = eg\ + e 3 ^ 2 <?2 + • • • for it, v±, b x , 6j_ and P. 

Substituting the power series expansions in the dissipative MHD equation we obtain 
in the first order approximation a solution that recovers the results of the linear theory. 
The compatibility condition for the equations of the second order approximation gives 
the governing equation for the wave motion in the slow dissipative layer, 

Ac{x-x c ) dv\\ di >\\ , d2 h c 2 Tc co S a dP 

2 ~m v — ~de +Av w ~d¥~ Xl aF-~ - ^ ^' (37) 



where 



Z\ c = -^-(l/ 2 -4 cos 2 a) , (38) 



dx 



A _ v Ac[(l + 1 ) v Ac + ^ c Sc\ cos a _ 770 c Tc r; 

\ v Ac + c sJ Wc "Ac 

When deriving Eq. (|37p we considered that vu » e 1 / 2 -D 1 1 2_ . In Eq. (|37p P is considered 
as a given function and it is determined by the solution to the linear ideal MHD 
equations outside the dissipative layer. 

As we have already mentioned, the main idea of connection formulae is to con- 
sider the dissipative layer as a surface of discontinuity and use the expressions for 
the jumps of P and u across the dissipative layer as the boundary conditions for 
the system of eq uations (I2 9 JI. The jump of total p r essure is the same as in linear 



theory (see, e.g., Hollwee Il988l : ISakurai et al.l Il99ll ; iGoossens and Rudermanl 1 19951 : 



iRuderman and Goosscnn 



[P] = 0, (40) 
where, in general, the jump of a function / across the dissipative layer is defined by 

[/]= lim {/(fl -/(_£)}. 
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(recall that £ = e 1 ' 2 (a; — i c )). To calculate the expression of [u] we use the relation 
between u and wn obtained in the first order approximation in the process of derivation 
of Eq. J37]), i.e. 

ox v* A 08 

where we have assumed u ~ eu\ and w e 1 ^ 2 i3|| x . It immediately follows from Eq. (|41[) 
that 

4-c cos a f°° dv\ 



u = - 



V 2 

v Ac 



V ^dx. (42) 



In the above equation V, denotes the Cauchy principal part of the integral which is 
used since v» ~ x~ l as x — > 00, so that the integral in Eq. (|41f) is divergent. 

In linear theory the jump in the normal component of velocity, [«], can be obtained 
explicitly in terms of P. In contrast, in the nonlinear theory we, in general, cannot 
solve Eq. (|37[) analytically and obtain an explicit expression for [u], which makes the 
connection formulae perhaps less practical. However we will see in Sect. [5] that there is 
on e exception. 

iBallai et al.l ( 1998ah extended the derivation bv lRuderman et al.l ( 1997af ) where the 



dominant dissipation processes are compressional viscosity and thermal conduction 

along the magnetic field lines. They formally assumed that N a ~ 1, so that R a ~ 

1/2 

e ' . Only the motion perio dic with respect to 8 has been considere d . Using the same 
procedure as one adopted by iRuderman et al.l ( 1997al ~). IBallai et al.l (1998a) obtained 



where now 



that the governing equation for the motion in slow dissipative layers is 
A c (x-x c ) dv\\ My 9 2 5|, 4 C cos a dP 

— v — W " W + Xa ~dW = 7^17 ~m ' (43) 

_ m V\2v\ c + 2cl L ) 2 (7-l) 2 m pK|| F 2 
Aa — 7. 9 — 9 — T~o 9~ r ' ; — 5 • K^l 

3Pc v Ac C Sc( v Ac + c Sc) 1P0ck B C Z Sc 

The main difference between Eq. Q44[l and Eq. (|37|) is in the terms describing the 
effect of dissipation, which is the last term on the right-hand side of Eq. (|44p . where 
it is proportional to the second derivative with respect to 8, while in Eq. (J37J) the 
corresponding term is proportional to the second derivative with respect to x. Another 
difference is that Eq. ()44|l does not con t ain the derivative with respect to time. This 
is related to the fact that IBallai et al. (1998a) did not allow slow time variation of 



the wave shape in the dissipative layer. The extension of the derivation to include this 
effect is straightforward. 

The derivation of governing equations for wav e motion in slow dissipati ve layers has 
been extend ed to include the e ffect of steady flow (|Ballai and Erdelyill 19981), cylindrical 



equil i brium (IBallai et al.l]200(]h and the twist of magnetic field lines ( Ballai and Erdelvll 
120021 : lErdelvi and BalIaiH2002 



Recent 1 v Tc4^ck~and" Ballail (|2008h generalized the derivation of governing equation 



for the wave motion in slow dissipative layers in order to take the Hall effect into 
account. The ratio of the second term on the right-hand side of the generalized Ohm's 
equation H14|l describing the Hall effect to the first term describing the plasma resistivity 
is of the order of io e r e . In the solar corona typically uj e T e > 10 6 , so that the Hall term 
strongly dominates over the resistive term. In classical MHD the magnetic field is frozen 
in the plasma when the plasma is infinitely conducting. When the Hall term is taken 
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into account the magnetic field is frozen in the electronic component of an infinitely 
conducting plasma, while the ions can drift with respect to the magnetic field lines. 
Hence, the account of the Hall term is equivalent to the account of ion inertia. The 
M HD theory where the Ha ll effect is taken into account is called the Hall MHD. 

IClack and Ballail (|200Sl ) repeated the derivation carried out bv lBallai et al. (1998a) 



however in the framework of the Hall MHD. As a result they obtained that the dynamics 
of waves inside the dissipative layer is given by 

Ac{x-xc) dv\\ 9v\\ , 5 2 S|| «3|| M|| 4 C cos a dP 

v W - Av " W + Xa ~bW + ha ~K W = 7^57 IS ' (45) 

where 

7i(u} e T e )k 2 V tana 

h - = -2 • ( 46 ) 

c Sc + v Ac 

In case when the dominant dissipative processes are the compressional viscosity and 
thermal conduction the characteristic thickness of slow dissipative layers given by the 
linear theory is 

x a kVX a ,.„.. 

Sc = w (47) 

This result is easily obtained assuming that the first and third terms in Eq. (|4]) are of 
the same order. Using Eq. (|47[) we obtain that the ratio of the second nonlinear term 
proportional to h a in Eq. Q to the first nonlinear term proportional to A is of the 
order of 

P{u eTe ) ( P_ + J_ 



Rm \ Re Re 

where the plasma-beta is defined by ft = Cg c /v\ c . In particular, in the corona, where 
P < 0.01, P e ~ 100 < R%, R m = 10 12 ~ 10 14 and cj e r e ~ 10 7 , this ratio does not 
exceed 10 -5 . Probably, it can be more pronounced in the chromosphere where aj e T e is 
already quite large, but R m is much smaller than in the corona. 



5 Explicit connection formulae 

As we have seen in the previous subsection, in general, the first connection formula, 
which expresses the continuity of the total pressure, is given in the explicit form no 
matter what is the value of the nonlinearity parameter in a slow dissipative layer (see 
Eq. (|40l) ). In contrast, we cannot write down an explicit expression for the jump of 
the normal component of velocity [u] , in general. Instead, to find [u], we first need 
to solve the governing equation for vu , then substitute Cm in Eq. ()42|l and calculate 
the integral in this equation. As a result, we cannot separate solutions in the resonant 
layer and in the external regions. Instead, we have to solve the linear ideal MHD 
equations in the external regions and the governing equation for in the dissipative 
layer simultaneously using the connection formulae Eqs. (|40[) and ()42[) to connect the 
solutions. This makes solving any problem involving nonlinear slow resonance quite 
complicated. 

However, there are two exceptions when we can obtain the second connection for- 
mula in an explicit form. The first one is the linear theory. In this case the con nection 
formulae for the cylindrical geometry have been derived bv lErdelvi et~aT1 (| 19951 ). These 
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formulae are easily translated to the planar geometry. We do not give these formulae 
here because our aim is to study nonlinear effects in dissipative layer, but refer to 
Goossens et al. (2010) in this volume. 

The second exceptional case is the limit of very strong nonlinearity, when non- 
linearity in a slow dissipative layer strongly dominates over dissipation ( Nj > 1 or 
N a S> 1) . In this case the explicit expression for [u] has been derived by IRudermanl 
(|2000l ). In what follow s we briefly outli ne the derivation of this expression and the 
main points. Note that IRudermanl ( 2000l ) assumed that the equilibrium magnetic field 



is in the z-direction. Here we consider the general case when the equilibrium magnetic 
field has both y and ^-component, the angle between the z-axis and the equilibrium 
magnetic field being a. 

In what follows we consider that the dominant dissipative processes are the isotropic 
viscosity and plasma resistivity. We assume that the motion is periodic with respect 
to the coordinate z. so that the wave motion in the dissipative layer is described 
by Eq. (|37p with the time derivative equal to zero. Let us introduce the modified 
nonlinearity parameter 

N = e 3/2 {kx ) 2 Ri, (48) 

where xq is the characteristic spatial scale of inhomogeneity, k — 2-n/L, and L is the 
period with respect to z. The condition that nonlinearity dominates over dissipation in 
the dissipative layer is written as N 3> 1. Let us introduce the dimensionless variables 

# = k a - i x ~ x c)sign{ A c ) n _ VAv\\ ^_ 2c^ c APcos 2 a 



5 c (27V)!/3 ' * c |4 c |(2JV) 1 /3' "* P0c v 2 Ac \A c \S c (2N)V3' 

(49) 



where 



(50) 

_k\A c \ Vpoc v 2 Ac 

is the characteristic thickness of slow dissipative layer given by the linear theory. In 
the new variables Eq. (|37[) with dvn/dt = is rewritten as 

dq dq 1 d 2 q dQ , . 

2a M +2q dV-N^ = -^- (51) 

Now q is a periodic function of $ with the period 27r. Since = it follows (q) — 0. 
It is easy to see that this restriction on q is compatible with Eq. (|5ip . Let us look for 
the solution to Eq. (|5ip satisfying this condition in the form of asymptotic expansion 
with respect to the small parameter 1/N, 

oo 

q= J2 N ~ n <l {n) - (52) 
n=0 

Let us assume that the function Q(t?) (and, consequently, -P($)) takes its maximum 
value in the interval [0, 2-n] at exactly one point "&m- We denote this value as Qm and 
the maximum value of P as Pm- This assumption implies that Q(~&m + 2- ku) = Qm, 
and Q ("&) < Qm when ■& ^ + 2mi, where n is any integer number. IRudermanl 
(|200d ) has shown that, under this condition, there is a quantity oo > so that the 
function q^ can be written as 

, n , f -a + [Qm - Q(#)] 1/2 , M < < «?s, 

q {0) = <^ (53) 

I -a - [Qm - QW] 1/2 , $s < # < ^ + 
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for | a | < cro- The quantity is a function of a implicitly defined by 

1 f&S i r-2ir+-d M 

a = D{$ 3 ) = — / [Q M - Qm 1 ' 2 dti-^r [Qm - Q(#)] 1/2 d&. (54) 

Equation H54|) determines c/ -' for i?^/ < ■& < i?^/ + 27T. We use the periodicity of c/ ' 
with repsect to ■& to define it beyond this interval. When |<j| > cro, the function g*- ) 
satisfies the condition 

g (0) (-a,#) = -q {0) {a,$). (55) 

iRudermanl (|2000l ) has shown that, in the zeroth order approximation, the wave 
energy dissipation in the resonant layer occurs in slow MHD shocks inside this layer 
rather than in the whole layer as it occurs in the linear theory. The shape of shocks is 
defined by the equation •& = $s(<j) + 2irn, where n is any integer number. The typical 
dependence of on ■& at fixed a, \o\ < cro, is shown in Fig. [T] Figure [2] displays the 
typical shape of slow shocks in the resonant layer. With the aid of Fig. it is easy to 





Fig. 2 The typical shape of slow shocks in the resonant layer (the figure was adapted from 
the study by Ruderman 2000). 
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Once again, this equation defines for < ■& < -d^j + 2ty, while q^is extended 
periodically beyond this interval. 

Now we are in a position to calculate [u]. Using Eq. (|49|l we rewrite Eq. (|42|) as 

r-i- ^1\a c \ { 2N)^ a v r°° q(0)da (57) 



v\ c A W 



We note here that, due to the relation Eq. (|55|l . the integrals of q^ over ( — 00, — cjq) 
and (00,00) cancel out each other. Then, using Eqs. (|54[) and (|56|l . we obtain 

/OO p(TQ 
q (0) da= / g (0) da = -2D(tf)[Q M -Q(tf)] 1/2 
-00 J —an 



-[Qm-QW] 1/2 \ / [QM-Qm 1 ' 2 ^ 



7i~ 



y [q m -qw] i/2 ^|. 



(58) 



Substituting this result into Eq. (|5Tf) and returning to the original dimensional vari- 
ables, we eventually arrive at 



2kV 3 c^ c d 
po c v\ c \A c \ dd 



[Pm-P(0)} 1/2 

2-7T+ 



J g [Pm - P(0)] 1/2 d9 - ^ [P M - P{6t' 2 (59) 

This is the second connection formula for a strongly nonlinear slow resonant layer. 
Recall that Pjy is the maximum value of P, Pm — P{9m)> an d k = 2n/L, where L is 
the period with respect to 8. 

Although we have used Eq. (|37[) to derive Eq. (|59[) . the derivation and expression 
for [it] are exactly the same in the case where the dominant dissipative processes are 
the compressional viscosity and thermal conduction, and the plasma motion in slow 
dissipative layers is described by Eq. (J43J) . 

On the basis of analysis presented in this section we can make one very important 

qualitative conclusion. The linear theory predicts that the dimensionless amplitude of 

1 /3 

large variables in slow dissipative layer is of the order of eR^ in plasmas where the 
dominant dissipative processes are isotropic viscosity and resistivity, and of the order 
of eR a in plasmas where the dominant dissipative processes are the compressional 
viscosity and thermal conduction (recall that e is the dimensionless amplitude of wave 
motion far from the resonant position). Hence, in accordance with linear theory, the 
amplitude of wave motion in slow dissipative layers tends to infinity in the limit of 
vanishing dissipation. 

However, it follows from the analysis in this section that, in accordance with non- 
linear theory just outlined, the amplitude of wave motion in slow dissipative layers is 
of the order of e ' in the limit of very small dissipation. Hence, nonlinearity causes 
saturation of the wave amplitude growth when the dissipative coefficients tend to zero. 
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6 Nonlinear effects in Alfven dissipative layers 

Resonant Alfven waves received much more attention than their slow wave counterparts 
since for a considerable time it was thought that Alfven resonance is much more im- 
portant in applications to the coronal low-beta plasma then slow re sonance. In fact the 
hunt f o r finding A lf ven w aves has not yet even been settled, see e.g. lErdelvi and Fedunl 
( 20071 ); Ijess et all (|2009l ). In the last few years the Alfven resonance received a new 



connotation related to the rapid damping of coronal loop kink oscillations, which is 
attributed to the resonant coupling of kink oscillations and local Alfven waves (see the 
review by Andries et al. 2009). 

Alfven waves are incompressible and transversal, therefore the dissipative mecha- 
nisms affecting Alfven waves are the shear viscosity and magnetic resistivity. Under 
coronal conditions the coefficients describing the magnitude of these mechanisms are 
very small, however dissipative terms in the momentum and induction equation can 
become as large as other terms in the region containing large spatial gradients. 

Since dissipation is only important in a narrow layer embracing the resonant surface, 
the dynamics of waves outside the dissipative layer is described by the same system of 
linear MHD Eq. Q29p as in the case of slow waves. 

The Alfven resonant position, xa, is defined by 

V — v a(xa) cos a. (60) 

In accordance with the linear theory the characteristic thickness of the dissipative layer 
embracing the ideal resonant position is given by 

where now V — va cos a, and 

(62) 



d /T ,2 2 2 n 

A = dx^ ~ VACOS a ' 



and the subscript A indicates that a quantity is calculated at x = xa- The estimate 
2. /3 

8a ~ LR A immediately follows from Eq. (|61[) . where 1/Ra = 1/Re + 1/Rm- In all 

applications to solar physics Ra is very large ranging from about 10 6 in the photosphere 

to up to 10 14 in the corona. 

Earlier in Sect. [3] we discussed that the large variables in Alfven dissipative layers 

— 1/3 

are vj_ and bj_ . The liner theory predicts that they are of the order of tR A , where e 
is the dimensionless amplitude. In the Alfven dissipative layer the ratio of the largest 

nonlinear terms in the dissipative MHD equations to the largest dissipative terms is of 

— 2/3 

the order of Na ~ € Ra • ® n ^ ne basis of this estimation we should conclude that 
nonlinearity starts to compete with dissipation as soon as Na ~ !• However, if we try 
to derive the governing equation for wave motion in Alfven dissipative layers using the 
same procedure as one adopted to derive Eq. (|37fl . we obtain a linear equation 

,V r ,A | V(Vi+ POAV) d 2 vj_ _ V sin a dP 

{ A) 09 + P0A Aa dx* - poaAa d6 ' [M) 

where vj_ is the oscillatory part of vj_. The fact that Eq. (|63[l does not contain nonlin- 
ear terms despite considering the full nonlinear MHD system of equations is in stark 
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contrast to the nonlinear description of wave motion in slow dissipative layers pre- 
sented in the previous sections. In fact. Iciack et al.l ( 2009aT l have shown that the linear 

1/3 

description of wave motion in Alfven dissipative layers remains valid if ei?^ -C 1, 
i.e. if the dimensionless amplitude of wave motion is much smaller than unity. Hence, 
the nonlinear description is only needed when the wave motion amplitude becomes of 
the order of unity. Once again note a sharp difference between the slow and Alfven 
dissipative layers. In slow dissipative layer the motion becomes strongly nonlinear as 
soon as its amplitude becomes of the order of e ' . 

The result that the linear description is valid for the wave motion in Alfven dis- 
sipative layers is the motion am plitude is much smaller than unity was predicted by 
Goossens and Rudermanllj 19951 ). aut no explanation for this behaviour has been given. 



Clac^ltn^I3ariail|2009aF showed that this result is related to the fact that the largest 



nonlinear terms in Alfven dissipative layers cancel out each other. Using the series 

— 1/3 

expansions with respect to R . they also showed that the main nonlinear effect in 
Alfven dissipative layers is generation of magnetosonic waves. The ratio of amplitudes 
of these magnetosonic waves to the amplitude of Alfvenic motion is of the order of 
eR 1 / 3 < 1. 

We do not describe the resu lts of linear theory for Alfv en dissipative layers. They 
can be found elsewhere (see e.g. iGoossens et al.lll995l . l20ld ). 



7 Resonant interaction of externally driven waves with inhomogeneous 
plasmas 

The process of resonant absorption involves the interaction of two oscillating systems 
which results in the energy transfer between them. In this context, probably, the most 
obvious case is the resonant interaction between an external laterally driven wave and 
the local oscillations of the plasma. Resonant interaction between global and local 
oscillations takes place if the frequency (or phase speed) of the global wave matches 
any value from the local slow or Alfven continuum, as stated in the Introduction. 

The effect of nonlinearity in a slow dissipative layer on the int eraction of sound 
wave w ith an inhomogene o us mag netized plasma was first studied bv lRuderman et al.l 
(19973) and iBallai et al. | |l998bh in the approximation of weak nonlinearity in the 



dissipative layer, and bv lRudermanl (|2000l l in the approximation of strong nonlinearity. 
In what follows we briefly outline the analysis carried out in these papers. 

We consider the equilibrium schematically shown in Fig. [3]described with use of the 
Cartesian coordinates x,y,z. The equilibrium magnetic field is unidirectional, parallel 
to the y2-plane, and makes the angle a with the z-axis. The inhomogeneous plasma 
layer (region II) occupies the slab < x < x^. It is sandwiched by two semi-infinite 
regions containing homogeneous plasmas (regions I and III), region I being magnetic- 
free and region III being penetrated by homogeneous magnetic field. Obviously the 
model we use here is very simplistic. In reality the magnetic structures are much more 
complicated. However the approximation we have made will help us to understand 
the fundamental characteristics of resonant absorption and its effectiveness. The equi- 
librium quantities in region I, II and III are labelled by subscripts 'e,' '0' and 
respectively. 
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Fig. 3 The sketch of the equilibrium state. The semi-infinite regions I (x < 0) and III (x > xq) 
contain homogeneous plasmas penetrated by homogeneous magnetic field, while the plasma 
and magnetic field in region II (0 < X < a;o) arc inhomogcneous. The resonant surface at 
x = x c is inside the dissipative layer shown by the shaded strip (adapted from Rudcrman et 
al. 1997b and Ballai et al. 1998b). 



All equilibrium quantities in region II depend on x only. The equilibrium quantities 
are assumed to be continuous and satisfying the condition of the total pressure balance, 

p e =p ( x ) + ^M=p i + ^L. (64) 

This equation, in particular, implies that the ratio of equilibrium densities in regions 
III and I is given by 

— = 2c ^ e (65) 
Pe 2c% i+1 v 2 M 

The plasma dynamics outside the dissipative layer is described the system of linear 
ideal MHD equations. In what follows we consider solutions in the form of propagating 
waves with permanent shape, and assume that perturbations of all quantities depend 
on x and 9 = z — Vt. Then the system of linear ideal MHD equations can be reduces 
to the system of two equations for the total pressure and x-component of the velocity, 

du = vop op _ Po D A du 

dx F 9(9 ' dx V 86' [ ' 

where the quantities F and Da are given by Eqs. (|30|l and (|31[) . This system is written 
for region II. To use it in regions I and III we have to substitute po by p e and pi 
respectively. 

The physical picture of the wave interaction with the inhomogcneous plasma is 
as follows. The sound wave incoming from region I interacts with the inhomogeneous 
plasma in region II. It is partially reflected back in region I, partially penetrate in 
region III, and also partially absorbed in the dissipative layer. As a result we have 
an outgoing (or reflected) wave in region I in addition to the incoming wave and a 
transmitted wave in region III. 
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In what follows we assume that the incoming sound wave is monochromatic. We 
write its wave vector as k = 0, k). Then its frequency is given by 

u? = fc 2 4 e (l + X 2 ). (67) 

The solution to Eq. (|66p describing the incoming wave can be written as 

P = e P eCOs[k(6 + xx)], u = cos[k(e + X x)}, (68) 

peV 

where V = uj/k = c<j e (l + X 2 ) 1 ^ 2 - In general, the outgoing wave will contain not 
only the fundamental harmonic, but also the overtones. Hence, the expression for the 
pressure perturbation in the outgoing wave can be written as P = ep e A(6 — X x ), where 
A is the function to be determined. Then the general solution to Eq. (I66|l describing 
the wave motion in region I is 

P = ep e {cos[k(6 + X x)\ + A(0 - X x)}, u = e^{cos[k(6 + X x)] - A(9 - X x)}. (69) 

pe v 

In region III the system of Eq. (]66f) can be reduced to 

2L< 29P a V 4 -V 2 (c 2 , l + v 2 Ai ) + v 2 Ai c 2 Sl cos 2 a 

g x +K ige U, S (c| i + «2. )(1/2 _ c 2. cos 2 Q) ■ 

In what follows we assume that there is a resonance position x c in region II, so that 
V — ct(x c ) cos a. In addition we consider that ct{x) is a monotonically increasing 
function in region II. Then it is straightforward to show that k 2 > 0, so that the wave 
motion in region III is evanescent. Using this result and expanding P in the Fourier 
series with respect to 9 we can write the solution to Eq. (|66[> decaying as x — ¥ oo in 
the form 

oo 

P = S~] P„exp[A;(wi0-«i|ra|(a:-xo)], (71) 



where, at present, P n are arbitrary constants. Substituting this result in Eq. (|66|l we 
obtain that the ^-component of the velocity in region III is given by 



Pi{V 2 - v 2 A , cos 2 a) 



P n sign(n) exp[k(in9 - nAnUx - x ))]. (72) 

. . cos z a ) L — ' 



To obtain the solution in region II we expand P and u in the Fourier series with 
respect to 8, and write the Fourier coefficients of these expansions as 

Pn(x) = p e V 2 X n (x), U n (x) = iVY n (x). (73) 

Substituting these expressions in Eq. (|66|l we obtain the system of ordinary differential 
equations for X n and Y n , 

^ = - ^ V2 - V } 2 COs2a K kY n , ^ = £?XlnkX n . (74) 
dx peV dx F 

Consider the two solutions to this set of equations satisfying the boundary conditions 
at x — 0, 

*i„ = l- Y 1 ~=0, (75a) 
*2n=0, Y 2 -=sign(n), (75b) 
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and the two solutions satisfying the boundary conditions at x = airj, 

XL = 1- Yin = 0, (76a) 
A 2 +=0, K+=sign(n). (76b) 

The solutions satisfying boundary conditions (|75p are the two linearly independent 
solutions to the set of equations (|74|l regular in the interval [0,x c ). Any other solution 
regular in this interval can be written as a liner combination of these two solutions. 
Similarly, the solutions satisfying boundary conditions (|76|l are the two linearly inde- 
pendent solutions regular in the interval [x c ,xq\. Any other solution regular in this 
interval can be written as a liner combination of these two solutions. 

Since the total pressure and rr-component of the velocity are continuous at x = 
and x = xq, the quantities X n and Y n satisfy the boundary conditions at x = 0, 

v - epe 

-K-n — 



(!<*1M + A «) . F « = (l 5 l\n\ - An) , (77) 

and at x = x c , 



_ Pnjxp) _ Ksign(n)P n (x Q ) 



Pe^ 2 ' ft(^2 



IT 



where <5jj is the Kronecker delta-symbol, and A n and P n are the coefficients in the 
expansions of functions A(ff) and P(x, 8) in the Fourier series with respect to 6. Using 
these boundary conditions we immediately obtain that 



An — 



Y„ 



PeV 2 



in x < x c , and 



+ A ") X ln -2Xsign(n) (§<*i|„| - A n ) X 2 „] , (79a) 
(j s i\ n \ + ^-*xsign(n) (j*i| n | --A»J Kjn] » ( 79b ) 

^A+), F n = ^(F+- ? y+), (80) 



y _ Pnjxp) f + + \ _ -Pn(ao) 



in x > x c , where 



«:p e U 2 



(81) 



pi{v 2 M COS 2 Q - V 2 ) 

It follows from the first connection formula, Eq. (|40[) . that X n (a;) is continuous at 
x = a^c- 

Equations (|79p and (|80p give the solution in region II. This solution contains the 
Fourier coefficients of the unknown function A(8). Now we are in a position to derive 
the equation for A(8). For this we calculate the jump of u across the dissipative layer 
using Eqs. ()79|) and (|80|) . and then compare it with the jump of u given by the second 
connection formula. As a result we obtain the integral equation for A(ff). 

Now we need to make one final remark. The procedure of derivation of the in- 
tegral equation for A{6) described above is only valid when the equilibrium mag- 
netic field is in the z-direction, i.e. when ot = 0. The reason is the following. We 
have assumed that or(x) is a monotonically increasing function and the slow reso- 
nant position is determined by the equation cj<(xc) cos a = V. Since va > or, we 
have that va(x c ) cos a > V. Since va(0) = 0, this implies that there is at leas one 
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point xa £ (0, x c ) where va(xa) cos a = V, i.e. there is at least one Alfven reso- 
nant position in (0, x c )- There is exactly one Alfven resonant position when va(x) is 
a monotonic function, while there could be a few Alfven resonant positions if va(x) 
is non-monotonic. In the simplest case when there is exactly one resonant position, 
xa < x c , we have to solve the ideal linear MHD equations in three intervals, (0,2:^), 
(xa,x c ) and (x c ,xq), use the solution to calculate the jumps of u at xa and x c , and 
then compare these jumps with those given by the connection formulae at the slow and 
Alfven resonance. The case where a = is exceptional because in this case there is no 
Alfven resonance in spite that there is xa satisfying va(xa) = V. 



7.1 Approximation of weak nonlinearity in dissipative layer 



Ruderman et all (jl997bh studied the interaction of sound wave with an inhomogeneous 



plasma where the dominant dissipative processes are resistivity and isotropic viscosity 
under the assumption that a = 0. In the nonlinear theory, in general, we cannot obtain 
an exp licit expression for [it]. To make analytical progress, however. iRuderman et all 
(|l997bl ) assumed that nonlinearity in the slow dissipative layer is weak and considered 
the nonlinear term in Eq. (|37|) as a perturbation. Then they carried out the asymptotic 
analysis using the modified nonlinearity parameter N given by Eq. (|48[) as a small 
parameter. In particular, they calculated the coefficient of resonant absorption defined 
as 

K= n in - gout ^ (g2) 

where 77; n and 77 ou t are the energy fluxes of the incoming and outgoing sound waves, 
respectively. They found that 

A « A'l + TV 2 A COT , (83) 

where Al is the coefficient of resonant absorption given by the linear theory, and 
A 2 A CO r is the nonlinear correction, and obtained the analytical expressions for Al 
and Acor in the approximation of thin inhomoge neous layer (fczp * C 1). W e do not 
give these expressions here. They can be found in IRuderman et all (|l997bl ). We only 
mention that A CO r < 0, so that nonlinearity supressess resonant absorption in the thin 
inh omogeneous l ayer ap proximation. 

iBallai et al.l ( 1998m carried out a similar analysis but for plasmas where the dom- 



inant dissipative processes are compressional viscosity and thermal conduction, once 
again assu ming that a = 0. The y used (kxo)N a as a small parameter, where N a is given 
by Eq. (t27)l . IBallai et al ] (|l998bh obtained the expression for A similar to Eq. ([83]), how- 



ever with (kxo)Na substituted for TV. Once again, A CO r < when kxo <C 1, however the 
ex pression for A CO r i n term s of equilibrium quantities i s different from that obt ained 
by Ruderman et al. | (|l997bh . The results obtained by IRuderman et aL (1997b) anc 



IBallai et all (|l998bh clearly show that, while Al remains the same no matter what 



type of dissipation operates in the slow dissipative layer, the coefficient of resonant 
absorption given by the nonlinear theory does depend on the dissipation type. 

In a more recent study Clack and Ballai (2009a) investigated the efficiency of 
resonant absorption when inside the dissipative layer the dispersion was of the same 
order of magnitude as nonlinearity and dissipation. Using the limit of weak nonlinearity 
these authors found that the effect of dispersion is to increases the absolute value of 
the nonlinear correction (see Eq. l83[) , i.e. to decreases the net coefficient of absorption. 
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Despite the change in the nonlinear correction this term still remained small compared 
to its linear counterpart. 



7.2 Approximation of strong nonlinearity in dissipative layer 

In this subsection we review studies about the same problem of absorption of the 
incoming wave energy in a slow resonant layer, however under assumption that the 
wave motion in the slow dissipative layer is strongly nonlinear. To simplify the analysis 
we assume that the equilibrium magnetic field is in the z-direction (a = 0), and region I 
is magnetic-free. The latter assumption implies that th e incoming wave i s a sound wave. 
The detailed investigation of this problem is given by IRudermanl J2000l). Here we onl y 



outline his analysis. In what follows we use the same notation as in iRudermanl (|2000h 
The solution to the linear ideal MHD equations in regions I and III, and in region II 
outside the slow dissipative layer, have been already described above. As it has been 
explained, in order to obtain the equation for function A(ff) describing the outgoing 
wave we have to calculate [it] using the solution of ideal linear MHD equations outside 
the dissipative layer, and then compare it with the expression for [u] given by the 
second connection formula. In the case of strong nonlinearity in the dissipative layer 
this expression is given by Eq. ()59[) . As a result we arrive at the integral equation 
determining A(0), 

Mf cos(k6) + X M{ sm(fc0) + M[A{0)] 

(>L-\-6m 1 

2k({S M - S(0)]* I / [S M - S(0)]* dd - I [S M - S(6)]z d0 I . (84) 



Here 

S(*) = ^f, S M = ^, C=^^, (85) 

ep e epe TTp 0c v\ c \A\ 

Recall that P^j is the maximum value of P in the dissipative layer. The operator 
M[A(0)] is given by 

00 

M[A(6)] = V - \M n A n exp(ink0) + M* x A* n exp{-ink0)} , (86) 
n=l 

where A n are the Fourier coefficients of function A(0), the asterisk indicates the com- 
plex conjugate quantity, and the coefficients M{ , M\ and M n are expressed in terms 
of the equilibrium quantities and the solutions (-^^j^i^) an d (X^„ , ) to the set 
of Eq . (|74p . We do not given these expressions here, they can be found in IRudermanl 
((2000). The function S(0) is expressed in terms of A(0) by 

5(6>) = Af lc cos(A;6i) + x^Tic am(k6) + C[A(0)}, (87) 

where the operator £[.A(0)] is given by 

00 

£[A{B)\ = J2 [( X inc + exp(mfc0) + (Af n c - i X X^ nc )A* n exp(-ink0)\ . 
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Once again, the coefficients X^ nc and ^2n C are expressed in terms of the equilibrium 
quantities and the solutions to th e set of Eq. (17411. And, once again, we do not given 
these expressions here referring to iRudermanl (|2000h instead. 

Since the operators and .M^^)] are expressed in terms of the Fourier 

coefficients of the function A(8), and these Fourier coefficients, in turn, are expressed 
in terms of integrals of the function A(9), the operators C[A(9)] and «M[v4(#)] are 
integral operators. Hence, Eq. (|84|1 is the integral equation for the function A(9). 

In spite that Eq. (|84|1 being a very complicated nonlinear equation, it has an ex- 
tremely simple solution of the form 

A(6) = acos(k8 + <p), (89) 

where the quantities a > and if are expressed in terms of the equilibrium quantities 
and the solutions to the set of Eq. (|74|l . We see that the outgoing wave is monochromatic 
in spite that the nonlinearity in the dissipative layer generates higher harmonics. The 
reason is that the amplitudes of these higher harmonics damp at distances of the order 
of thickness of the dissipative layer, so that only the fundamental harmonic survives 
at large distances from the slow resonan t position. 

Using the solution given by Eq. flU) IRudermanl (|2000h calculate the coefficient of 



resonant absorption. It immediately follows from Eq. (|82[) that K — 1 — a . Note that 
this expression is valid regardless whether we use the linear or nonlinear description of 
plasma motion in the slow dissipative layer. The nonlinearity only affects the value of 
a. In general, the set of Eq. (|74[> can be solved only numerically. However, in the thin 
inhomogeneous layer approximation (kxQ <C 1) the analytical solution can be obtained 
in a straightforward way. In that case the expression for the coefficient of resonant 
absorption is given by 

^nl = 4^2 +O{k 2 x 2 ), (90) 

where we have used the subscript 'NL' to indicate that the coefficient of resonant 
absorption has been calculated using strongly nonlinear description of the wave motion 
in the slow dissipative layer. In Eq. (|90|l the quantity ? is given by Eq. (|81[) with a = 0, 
and £ is given by Eq. (|85[) . Note that ( = O(kxo), so that -Knl — O(kxo), which implies 
that resonant absorption is weak in the thin inhomogeneous layer approximation. For 
the ratio of the coefficients of resonant absorption calculated using linear and strongly 
nonlinear descriptions we obtain 

^ = 4 r +0(fcz o )«O.81. (91) 

We see that, in the thin inhomogeneous layer approximation, nonlinearity reduces the 
efficiency of resonant absorption. This result is in good agreement with the results 
obtained in the weak nonlinearity approximation and described in the previous sub- 
sec tion. 

IRudermanl ( 2000h studied numerically the dependence of -Knl/^L on kxQ for dif- 
ferent values of x an< A a typical equilibrium. He found that A'nl/^'l is a non-monotonic 
function of kx$. First it grows, takes its maximum at k m XQ, and then monotonically 
decreases. Typically k m XQ is between 4 and 6, and the maximum value of -Knl/^l is 
larger than 1. However the most important result is that, for < kxo < 10, -K^nl/^L 
does not deviate from unity by mode than 20%. Hence, we conclude that the effect 
of nonlinearity in slow dissipative layers on the coefficient of resonant absorption is 
moderate. 
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As pointed out previously, the nonlinear coefficient of resonant absorption depends 
on particular dissipative processes operating in a slow resonant layer. However, this 
dependence disappears in the approximation of strong nonlinearity. This result is re- 
lated to the fact that, in the approximation of strong nonlinearity, dissipation in a slow 
resonant layer occurs not in the whole volume of this layer, but in slow shocks inside 
the layer. It is a very well known property of shocks that, while their internal struc- 
ture is determined by particular dissipative processes operating in a shock, the amount 
of energy dissipated at the shock remains the same regardless what the dissipative 
processes are. 

Recently Clack et al. (2010) introduced the concept of coupled resonances within 
the framework of nonlinear resonant MHD, where the transmitted waves at one of the 
resonances could play the role of the incoming wave for the second resonance. It is 
obvious that due to the f3 <C 1 in the solar corona a coupled resonance (slow+Alfven) 
would be impossible, however in plasmas where the plasma-/? is much closer to unity 
the coupled resonance could take place. In order to ensure that the second resonance 
takes place before the transmitted waves becomes evanescent, the proximity of the two 
resonances must be smaller than k~ . When looking at the efficiency of the coupled 
resonances Clack et al. (2010) found that the absorption coefficient of the coupled 
resonance was larger than the sum of individual coefficients taken separately at the 
two resonances. 



8 Generation of mean flows 



When a wave propagates through a medium, the nonlinear interaction of the wave 
and medium generates a flow. This process is very well known in nonlinear acoustics 
where flows generated by pr opagating sound waves are called "acoustic flows" (see e.g. 
iRudenko and SoluvarJ[l977l ). 

We have already mentioned about the generation of flows by resonant MHD waves 
in Sect. [4] when we split the velocity in the slow dissipative layer in the mean and 
osci llatory parts . The f lo ws very importan t for res onant and diagnostic purposes (see 
e.g. iDovle et al] (|l997t )V iRuderman et all |l997aT ) derived the formulae determining 
the jumps of the derivative of the the mean flow across the dissipative layer under the 
assumptions that the plasma motion is periodic with the period L, and the dominant 
dissipative processes are isotropic viscosity and resistivity. Their study was later ex- 
tended by Ballai et al. (2000b) for cylindrical geometries. These formulae determining 
the jumps can be written as 



dv 



d.i: 



LV 



1 + 



P0cC Tc y 



/>/:(s 



dx, 



(92) 



dim 



dx 



(VO + POcVW 



V0V 2 Al 



L cos a 



dx. 



(93) 



When the motion in a slow dissipative layer can be described by the linear dissipative 
MHD equations (which is possible if N$ <C 1) we can find the explicit expression for 
the integrals in Eqs. (|92[1 and (|93|l in terms of perturbation of total pressure P. In that 
case «|| in the slow dissipative layer is described by Eq. (|37[) with the term proportional 
the time derivative and the nonlinear terms (which are the first and third term on the 
left-hand side) are neglected. This linear equation can be easily solved either by the 
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metho d used bv lGoossens et al ] (|l995l l or by the method used bv lTirrv and Goossens I 
(| 19961 ) to obtain 



ic Tc V 2 P n 
P0cV 2 Ac A c 



2irnA c 



XiLV 



1/3 



sign(n A C )F I (x 



2-kuAc 



■ x c 



XiLV 



1/3 



sign(n/A c 



(94) 



where £>||„ and P n are the coefficients in the expansions of tiii and P in the Fourier 
series with respect to 6, and the _F-function is defined by 



exp(ij/CT — a' /3) da. 



Changing the order of integration and using the Parseval identity yields 



POO 


9v\\n 


J — oo 


dx 



dx. 



(95) 



(96) 



Using Eqs. (|94l) and (|95[) it is not difficult to obtain that 



<:)b 



\\n 



2iY 2 c 2 Tc V i n\P n 



dx 



dx = 



P0c v Ac 



A C \LX, 



Substituting Eqs. (j96j and (J9TJ) in Eqs. ((92j) and (j93j we arrive at 



dx 



.22x7-2- 00 

477 C ^ V SmQ ^n|P„ 



dun 



4^0^^(770 + p 0c ?7) cos a / r? 



POc 



+ 



c TcV 
v 2 



-1 00 



(97) 



(98) 



(99) 



It is instructive to estimate the order of magnitude of these jumps. It is not difficult 
to show that both jumps are of the order of [V / L)Nf R e -R^ ^ . It is interesting that, 
although the amplitude of the wave motion in the dissipative layers is small, the jumps 
given by Eqs. (|98|) and f|99f) can be quite large if R e 3> Rm, i.e. if the magnetic Prandtl 
number, P m = Vo/POcVi ' s small. Simple estimates show that these jumps are larger or 
of the order of V/L when 1 R m <^ R e and R e > e~ 2 . 

To find the profiles of the components of the mean velocity we need to impose 
boundary conditions far away from the dissipative layer. For example, if there are 
rigid walls at x = ±a where the condition of adhesion has to be satisfied, then the 
components of the mean velocity take the simple form 



dv 1 



dx 
dvj_ 



2 ' 
x + a 



x < 0, 
x > 0, 



dx 
d«|| 
dx 



x < 0, 



x > 0. 



(100) 



It is worth mentioning that the derivations carried in cylindrical geometry by Ballai 
et al. (2000b) resulted in some analogue relations for the jump in the mean flow com- 
ponents, but the properties of the mean flow generated by the nonlinear resonant slow 
waves are the same. This means that the results presented here are rather robust. 
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Although the dynamics of waves at the Alfven resonance can be described within 
the framework of linear MHD with great accuracy, the process of resonance can produce 
a mean flow even in this later case. The problem of mean flow generation at the Alfven 
resonance was recently studied by Clack and Ballai (2009b). In this case the jump in 
the oscillatory part of he parallel and perpendicular components of velocity are given 
as 

,2-3 00 

Z^n\P n \, (101) 



dx 



and 



dx 



k 



2 sin 2 a cos a 



2up^\Aa\ 



(102) 



For the particular case of a = tt/4 Clack and Ballai (2009b) found that these jumps 
scale as e^'^V/L. For example, if the dimensionless amplitude e = O(10 -4 ) then the 
predicted mean shear flow is of the order of 10 km s _1 in both upper chromosphere 
and solar corona. This value is comparable with observed values of bulk flow in coronal 
loops. However, we should keep in mind that the important difference between these 
two types of flows is their range; mean shear flows are rather localised while bulk 
motion of plasma occurs often along the entire coronal loop length. 

Studies have been carried out to investigate the properties of shear flows, however 
nearly all of these have been numerical due to analytical complications when consid- 
ering nonlinearity, turbulence and resonant absorption simultaneously. These studies 
have found that shear flows c ould give rise to a Ke fvin-Helmholtz instability at the nar- 
row dissipative layer (see, e.g. lTerradas et al ] l2008h . This instability can drive turbulent 
motions and, in turn, loca lly enhance transport coefficients which can alter the effi- 
cienc y of heating (see, e.g. iKarpen et al.lfl994l ; lOfman at al.lll994l ; lOfman and Davilal 
Il995l l. The generation of mean shear flow can supply additional shear enhancing tur- 
bulent motions. 



9 Conclusions 

Resonant absorption is a mechanism that ensures the effective transfer of energy be- 
tween interacting systems. In solar and space plasmas resonant absorption is used in 
conjunction with the interaction of global waves and oscillations with inhomogeneous 
plasmas. Although the initial purpose of resonant absorption to be the phenomenon ex- 
plaining the heating of solar corona seems these days a bit too optimistic, the advances 
of observational evidences of the last couple of years showed that resonant absorption 
is still a remarkable and brilliant mechanism on its own, able to explain a series of 
delicate effects occurring in the solar atmosphere and beyond. 

The present review highlighted the progress made in the field of nonlinear reso- 
nance in the last 15 years. Simple estimations presented in this paper show that near 
resonances the amplitudes of oscillations can grow considerably and only a nonlinear 
description is adequate to describe accurately the physics near resonances. Luckily, as 
shown in the present review, these nonlinear additives to the effectiveness of absorption 
(compared to the linear counterpart) is not so significant, in general nonlinearity tends 
to decrease the net coefficient of absorption. In the case of resonant Alfven waves, re- 
sults show that the linear approach can give rather accurate results, with no need of 
using a mathematically more cumbersome nonlinear description. 
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We also showed that nonlinearity can have additional effects to the change in the ab- 
sorption coefficient. Due to the nonlinear absorption of wave momentum in the vicinity 
of resonance, a shear mean flow is generated that is continuous across the layer con- 
taining the resonant surface, while its derivative has a jump. Simple estimations of the 
magnitude of the mean shear flow show that in the solar corona and upper chromo- 
sphere this mean flow is of the order of 10 km s _1 , a flow comparable with the existing 
bulk motion of the plasma in coronal loops. The generated shear flows can generate 
instabilities that can enhance locally the transport coefficients, i.e. the effectiveness of 
resonance. 

The analysis on the nonlinear resonant absorption reviewed by our paper contains 
several simplifying assumptions that helped us progress in the analytical study of the 
process of nonlinear resonance, ft is obvious that considerable advances in this field 
can be made only using numerical simulations. One way to extend the existing theories 
is the study of nonlinear absorption in 2-D geometries. Further, all presented theories 
supposed (and that is true for all studies of resonant absorption, including linear cases) 
that the interacting waves are monochromatic, however it is more likely that waves in 
nature do not depend on one single wavenumber, but they depend on a spectrum of 
values meaning that resonant absorption would occur for all those wavenumbers for 
which the resonant condition is satisfied (some encouraging numerical studies 
were carried earlier by Ofman and Davila 1996 and Ofman et al. 1998). 
It remains to be seen what subtle effects the mean shear flow has on the resonance, 
probably through numerical simulations. It also known through numerical inves- 
tigations (see e.g. Ofman et al. 1998) that the resonant absorption leads to 
the modification of the loop density structure in the nonlinear regime, thus 
affecting the location, dynamics, and evolution of the resonant absorption 
layer, as well as the heating of the loop. This aspect still needs further 
investigation for the cases presented in our paper. 

The abundance of high resolution observations and the extended numerical possi- 
bilities available will play an essential role in the progress of resonant absorption's study 
and its applicability in explaining other phenomena occurring in solar and space plas- 
mas predicting a bright future for this simple, yet fundamental effect in inhomogeneous 
space and laboratory plasmas. 
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